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We present a combined x-ray diffraction and infrared spectroscopy study on the phase behav- 
ior and molecular dynamics of n-hexadecanol in its bulk state and confined in an array of aligned 
nanochannels of 8 nm diameter in mesoporous silicon. Under confinement the transition temper- 
atures between the liquid, the rotator R/j and the crystalline C phase are lowered by approxi- 
mately 20 K. While bulk n-hexadecanol exhibits at low temperatures a poly crystalline mixture of 
orthorhombic f3- and monoclinic 7- forms, geometrical confinement favors the more simple ^^-form: 
only crystallites are formed, where the chain axis are parallel to the layer normal. However, the 
7-form, in which the chain axis are tilted with respect to the layer normal, is entirely suppressed. 
The /3- crystallites form bi-layers, that are not randomly orientated in the pores. The molecules are 
arranged with their long axis perpendicular to the long channel axis. With regard to the molecular 
dynamics, we were able to show that confinement does not affect the inner-molecular dynamics of 
the CH2 scissor vibration and to evaluate the inter- molecular force constants in the C phase. 

PACS numbers: 64.70. Nd, 61.46.Hk, 81.07.-b, 61.43.Gt 



I. INTRODUCTION 

The physical properties of condensed matter spatially 
confined in pores or channels of a few nanometer in di- 
ameter can differ markedly from the behavior in the bulk 
state. In particular, phase transitions can be entirely 
suppressed or significantly altered in comparison to their 
bulk counterparts^'^^^^^. Also the dynamics of condensed 
matter confined in mesopores, most prominently in the 
vicinity of glass transitions^i^iii^i^ii^iiiii^ii^ii^, can be af- 
fected markedly. 

Intimately related to these changes in the phase transi- 
tion phenomenology the architectural principles of molec- 
ular solids can substantially differ in the spatially con- 
fined state from the bulk state. This depends, however, 
sensitively on the complexity of the building blocks. For 
simple van-der-Waals systems, such as Ar and N2, a 
remarkable robustness of the bulk structures has been 
found for the solid state under confinemen t -"^^i^^ ^^. By 
contrast, the structural properties of pore fillings built 
out of more complex building blocks, such as linear 
hydrocarbons^^^^^'^^^^^'^^ or liquid crystals^^^^^ are very 
susceptible to confinement on the meso- and nanoscale. 
For example, a quenching of the lamellar ordering of 
molecular crystals of n-alkanes has been observed in tor- 
tuous silica mesopores of Vycor^^. However, in tubu- 
lar channels of mesoporous silicon this building principle 
of hydrocarbon molecular crystals survives, albeit a pe- 
culiar texture has been observed for the pore confined 
solids^: The long axes of the molecules and thus the 
stacking direction of the lamellae are oriented perpendic- 
ular to the long axis of the pores. 

Here we present an experimental study on a medium- 
length, linear alcohol C16H33OH, a representative of the 
1-alcohol series, imbibed in mesoporous silicon. We ex- 
plore the phase behavior of the confined alcohol by a com- 



bination of x-ray diffraction and infrared spectroscopy 
measurements. As we shall demonstrate, we profit in 
those experiments both from the parallel alignment of 
the silicon channels and from the transparency of the sil- 
icon host in the infrared region. 



II. EXPERIMENTAL 

The porous silicon samples used in this study were pre- 
pared by electrochemical etching of a heavily p-doped 
(100) silicon wafer ^ with a current density of 13 
in a solution composed of HF, ethanol and H2O (1:3:1 
per volume These conditions led to a parallel ar- 
rangement of non-interconnected channels oriented with 
their long axes along the <100> crystallographic direc- 
tion of silicon, which coincides with the normal of the 
wafer surface. After the porous layer had reached the 
desired thickness of 70 microns, the anodization current 
was increased by a factor of ten with the result that the 
porous layer was released from the bulk wafer under- 
neath. Using nitrogen sorption isotherms at T = 77 K, 
we determined a porosity of 60% and a mean channel di- 
ameter of 8 nm. The single crystalline character of the 
matrix was checked by x-ray diffraction. Transmission 
electron micrographs of channel cross sections indicate 
polygonal, rough channel perimeters rather than circu- 
lar, smooth circumferences^^. 

The matrix both for the infrared spectroscopy and 
the x-ray measurements were filled completely via 
capillary action (spontaneous imbibition) with liquefied 
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CieHaaOH^^. Bulk excess material at the surface was 
removed by paper tissues. 

Infrared spectra in a range of wavenumbers V from 
4000 to 800 cm~^ with a resolution of 1 cm~^ were 
measured with a Fourier Transform Spectrometer (FTIR 
Perkin Elmer System 2000). This range corresponds to 
frequencies from 3-10^^ Hz to 1.2-10^^ Hz (wavelengths 
from 10 fivn to 2.5 fim). For both the bulk material 
and the filled porous samples the same sample holder 
was used, i. e. a copper cell with two transparent KBr 
windows. In the confinement experiments the long 
channel axes were oriented parallel to the beam axis, i. e. 
perpendicular to the electric field vector. The sample 
holder was placed into a cryostat (a closed cycle refriger- 
ator CTI cryogenics. Model 22) allowing us to vary the 
temperature from 50 to 340 K. The temperature was 
controlled with a LakeShore 340 temperature controller 
with a precision of ±0.25 K. All IR-spectra that we show 
in the following were measured during cooling (typical 
cooling rates were of the order of 0.5 K/min). 
Heating scans show the same behavior except for the 
transition temperatures, which are some degrees higher 
(see below). 

For the x-ray measurements the sample was mounted 
on a frame in a sample cell consisting of a Peltier cooled 
base plate and a Be cap. The cell was filled with He gas 
for better thermal contact. The Be cap sits in a vacuum 
chamber, the outer jacket of which has Mylar windows al- 
lowing the passage of the x-rays over a wide range of scat- 
tering angles within the scattering plane (see Fig. [1]). 
But the set-up allowed practically no tilt with respect to 
the scattering plane. The temperature was controlled by 
a LakeShore 330 over an accessible range from 245 K up 
to 370 K. The measurements were carried out on a two- 
circle x-ray diffractometer with graphite monochroma- 
tized CuKq, radiation emanating from a rotating anode. 
The porous sheet was mounted perpendicular to the scat- 
tering plane. The two angles that could be varied were 
the detector angle 20 and the rotation angle uo about the 
normal of the scattering plane. The samples were stud- 
ied as a function of temperature by performing several 
<l>-scans. In this paper we concentrate on radial 26-uj- 
Scans in refiection geometry, i.e. along qp with <l>=0°, 
and in transmission geometry, i.e. along qg with <l>=90° 
(see Fig. [I]). 



A. Structure of bulk n-hexadecanol 

n-Hexadecanol, C16H33OH, is an almost rod-like 
molecule with a length of 22 A and a width of 4 A. The 
C-atoms of the backbone are in an all-trans-configuration 
so that they are located in a plane^^. 

At low temperatures n-alcohols form bi-layered crys- 
tals in two possible modifications: the so-called 7-form, 
i. e. a monoclinic structure as sketched in Fig. [2] (C|^ — 




FIG. 1: (Color online) The angular variable ^ — u — O with 
0° for scans along qp and $ = 90° for scans along qs) is 
indicated. Also shown are the incoming and outgoing x-ray 
beam for scans along qp (solid lines: wave vector component 
parallel to the pore axis, this means reflection geometry), and 
along qs (dashed lines: wave vector component perpendicular 
to the pore axis; this is transmission geometry). For scans 
along qp, the detector angle 20 and the rotation angle uo are 
represented. 



A2/a^^'^^), or the so-called /^-form, i. e. an orthorhom- 
bic structure as sketched in Fig. [3?^. In the 7-form, 
the molecules include an angle of 122° with the layer 
plane. Within the layers, they are close-packed in a 
quasi-hexagonal 2D array, described by the rectangular 
in-plane lattice parameters a and h (according to Ref*^ 
a = 7.42 A and b = 4.93 A holds, so that a/b = 1.5). 
There are two different alternating orientations for the C- 
C-plane of the backbone leading to a herringbone struc- 
ture (see Fig. [2]3). The /3-form exhibits an identical ori- 
entational order of the backbone, but the molecules' axes 
remain perpendicular to the layers as sketched in Fig.[3f^. 
In addition, gauche- and trans-conformation of the CO- 
bond alternate with molecules in this phase, while they 
are in an all-trans configuration in the 7-form. In gen- 
eral, the 7-form dominates at low temperatures for the 
even alcohols, while the /3-form is more frequent in odd 
n-alcohols^^^^^. For n-hexadecanol both the orthorhom- 
bic /3-form^^ and the monoclinic 7-form^^^^^ are reported. 
Depending on the preparation conditions it is possible to 
obtain a polycristalline mixture of the monoclinic 7- and 
the orthorhombic /3-form^^. 

Upon heating, the crystalline phase undergoes a tran- 
sition into a so-called Rotator- (Il)-phase i?//, which is 
schematically depicted in Fig. [H ^. This phase has a 
hexagonal in-plane arrangement with the c-direction per- 
pendicular to the cell base. The hexagonal arrangement 



^ For several alkanes, there also exists a Rotator- (I)-Phase Rj, 
where the molecules switch between two equal positions. 
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FIG. 2: 7-form of the crystalline low temperature phase of 
bulk C16H33OH (T < 310 K). The structure is monoclinic. 
The left sketch shows the orientation of the molecules with 
respect to the layer normal n, the right sketch the in-plane 
arrangement, i. e. a projection of the backbones into the a-b- 
plane. Compare with the /3-form sketched in Fig. [3] 




FIG. 3: /3-form of the crystalline low temperature phase of 
C16H33OH. In contrast to the 7-form (see Fig. [2^), the long 
chain axes are not tilted but parallel to the layer normal n, i. e. 
the structure is orthorhombic. Bulk C16H33OH can exhibit 
a polycrystalline mixture of 7- and /3-form (see Sec. Ill A|) . 
Confinement into nanopores leads to the /3-form (see below, 
Sec. lnTAl) . 



can be indexed with an orthorhombic cell with a ratio 
of rectangular basal lattice parameters of a/h = ^/3 
On a microscopic level the change in the center of mass 
lattice from the low-temperature crystalline phase to the 
rotator phase can be attributed to the onset of hindered 
rotations of the molecules about their long axes between 
six equivalent positions (the stars in Fig. [4)3). Further 
heating above 322 K leads to the liquid stated. 



III. RESULTS 

A. Structure of confined n-hexadecanol 

We have determined structures, phase sequences and 
transition temperatures of n-CieHaaOH confined in 
mesoporous silicon by x-ray diffractometry. The upper 
panel in Fig. [5] shows diffraction patterns along qp at se- 




FIG. 4: Structure of bulk n-hexadecanol in the Rotator- (II)- 
phase (for 310 <T< 322 K), a hexagonal arrangement. The 
right picture shows the perfect hexagonal lattice in the a-b- 
plane. Confined C16H33OH exhibits the same structure in 
its rotator phase, but in a different temperature range (see 
below. Table ini- 



lected temperatures while cooling. The appearance of 
a broad Bragg peak at 20 21° indicates solidifaction. 
Its position is compatible with the leading hexagonal in- 
plane reflection of the Ru phase. Upon further cooling 
a second peak at 20 24° shows up. This change in the 
diffraction pattern indicates an uniaxial deformation of 
the hexagonal lattice. Both reflections can be mapped 
on a 2D rectangular mesh characteristic of an uniaxially 
deformed hexagonal cell. The overall resulting pattern 
is, however, incompatible with the monoclinic structure 
of the low temperature bulk crystalline phase. 

Additionally to the qp-scans, we performed also scans 
for a variety of additional orientations of the scattering 
vector with regard to the long axis of the channels. These 
patterns differ markedly, which is indicative of a strong 
texture of the pore confined cystallized alcohol. It is no 
powder in the crystallographic sense. In particular, there 
are strong in-plane reflections and no layering reflections 
for scans along qp, while the qg-scans for the same sample 
show at least very weak reflections characteristic of a bi- 
layer stacking and only very weak leading in-plane reflec- 
tions (see Fig.[5|). An analysis of the width of the layering 
reflections yields a coherence length of 7 (±1.5) nm. 

As discussed in more detail in Refs.^^ and^^, the overall 
picture which emerges from these results can be summed 
up as follows: the alcohol molecules form orthorhom- 
bic structures with a bilayer-stacking direction along the 
c-direction. Within the bilayers (the a-b-plane), the 
molecules' backbones are untilted with regard to the 
stacking direction and the backbones are orientationally 
either fully ordered (in a herringbone fashion) or par- 
tially ordered, as known from the R/ phase of n-alkanes. 
The superlattice reflection characteristic of the full, her- 
ringbone type orientational ordering has been searched 
for and could weakly be detected at low temperatures. 
The degree of uniaxial deformation of the hexagonal cen- 
ter of mass cell, quantified by the deviation of the ratio 
a/b from its value in the hexagonal phase (>/3), also in- 
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dicates a full orientational ordered state (see Table I,^^). 
Thus, the diffraction data are compatible with the bulk /3 
modification discussed above. This conclusion is also sup- 
ported by an analysis of the infrared spectroscopy data 
sets presented below. 

More importantly, the peculiar dependency of the 
diffraction patterns on the orientation of the q-vector 
with regard to the silicon host indicate that the bi-layer 
stacking direction is perpendicular to the long axis of 
the channels and, consequently, that the long axis of the 
molecules is oriented perpendicular to the long axis of the 
channels (see Fig. [6]). At first glance, this finding may 
appear somewhat counter-intuitive. Albeit it can be un- 
derstood as resulting from the crystallization process in 
a strongly anistropic, capillary-like confined liquid^^^^-. 
It is a well established principle in single crystal growth 
that in narrow capillaries the fastest growing crystalliza- 
tion direction prevails over other directions and propa- 
gates along the long axes of capillaries^^. For layered 
molecular crystals of rod-like building blocks this direc- 
tion is an in-plane direction, which is perpendicular to 
the long axis of the rods. If this direction is aligned par- 
allel to the silicon nanochannels due to the crystallization 
process, it dictates a perpendicular arrangement of the 
molecules' long axes with regard to the long channel axis, 
in agreement with the diffraction results presented here. 
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TABLE I: Lattice parameters a and b of bulk and confined 
C16H33OH and the diagonal d = va^-F^ of the subcell (see 
Fig. I14[) . The confined data result from our x-ray measure- 
ments and the bulk data are taken from the literature^-. 
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FIG. 5: (Color online) Upper panel: Diffraction pattern of 
qp-scans of pore confined C16H33OH while cooling at selected 
temperatures. Lower panel: Diffraction pattern of qs-scans 
at selected temperatures while cooling. The leading in-plane 
reflections are barely visible. Moreover, at lower angles {26 < 
10) there are very weak and diffuse layering reflections. 



The temperature dependent diffraction study allows us 
to gain additional information on the relative stability of 
the different nanochannel confined phases. In Table [TTl 
we display the phase transition temperatures of confined 
C16H33OH as inferred from the appearance or disappaer- 
ance of characteristic Bragg peaks. There is a hysteresis 
between heating and cooling for both the fluid-R//- and 
the R//-C-transition (8 K and 3 K, respectively). Com- 
pared to the bulk data (see also Tab. the transition 
temperatures of pore confined C16H33OH are lowered. 
On cooling, the lowering is of the order of AT= 18 K for 
the fluid-i?//-transition and AT= 26 K for the Ru-C- 
transition. This observation is analogous to phase tran- 
sitions shifts reported for other pore condensates^i^^ . 

Furthermore, the temperature range of the confined 
Rjj phase, 14 K upon cooling and 19 K upon heating, is 
larger than that of the bulk material (12 K). Obviously, 
confinement stabilizes the orientational disordered Ru 



phase, similarly as has been found for n-alkanes^^ and 
for other orientational disordered, plastic phases under 
spatial confinement^. 

Since the pores were completely filled at higher tem- 
peratures, when hexadecanol is in its liquid state, the 
pore filling at low temperatures does not consist only 
of bi-layer crystals: the change of volume at the R//-C 
phase transition is about 10% (see Fig. [7j), so that there 
are voids and/or molecules that are not part of a bi-layer 
crystal. However, our experiments do not give us infor- 
mation about their spatial arrangement. 
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TABLE II: Transition temperatures as determined by x-ray 
experiments of confined and bulk CieHssOH^^. They agree 
with those inferred from IR - measurements (CH2-scissoring 
vibration), see Figs. [TOlandfTTI 
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FIG. 6: Sketch of the orientation of pore condensed 
C16H33OH in the porous sheet. Left pore: view on the rect- 
angular in-plane arrangement of the molecules. Right pore: 
bi-layered crystal; the long molecule axis are oriented perpen- 
dicular to the pore axes. At maximum two bi-layers fit into a 
pore. 



B. Molecular dynamics 

The dynamics of bulk-CieHaaOH has already been 
investigated in IR- measurements in the pasti^i^?. In 
order to show later on how the molecular dynamics 
is affected by spatial confinement on the nm-scale, 
we display some of our bulk spectra in the following. 
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FIG. 7: (color online). Temperature dependence of the area 
per molecule A for C16H33OH confined in porous silicon. 



Here we focus on two characteristic vibrations, the 
OH-stretching and the CH2-scissoring vibration. 

Figs. [8^) and [9^) show the bulk spectra of the OH- 
stretching-band in the respective phases (compare with 
Figs. [IH]). In the liquid state (above 322 K) the peak 
maximum is located at about 3345 cm~^. A decrease 
of temperature below 321 K yields a shift of the peak 
position to about 3325 cm~^ indicating the molecular 
rearrangement in the Ru phase. A further decrease of 
temperature below 310 K results in a splitting into two 
peaks at approximately 3310 cm~^ and 3220 cm~^. Con- 
fined C16H33OH shows a different behavior. There is only 
one peak in the whole temperature range, the position 
of which changes reflecting the transition between liquid 
phase and Ru phase as well as between Ru phase and 
C phase (see Figs. [HJd and [9)3). 

The fact that the OH-band of bulk C16H33OH splits 
at low temperatures while no splitting is observed un- 
der confinement confirms the structural differences al- 
ready observed in the x-ray experiment. For example, 
Tasumi et. al have studied bulk alcohols CnH2n+iOH 
from n = 11 — 37 using infrared spectroscopy^^, Ventola 
et al. alcohols with n = 17 — 20^^. Those alcohols 
showing at low temperatures (C phase) the monoclinic 
7-form, such as C16H33OH, exhibit the splitting of the 
OH-band, while those that take the orthorhombic /3-form 
show a single peak. This is due to differences in the 
spatial arrangement of the hydrogen bonds as well as in 
the distances of neighboring 0-atoms: in the crystalline 
7-form, where the molecule axis are tilted (see Fig. [2]), 
the molecules show an all trans conformation, and the 
intra-layer 0-distance {c:^ 2.74 A) differs from the inter- 
layer 0-distance (:^ 2.69 A). However, in the orthorhom- 
bic /3-form (Fig.[3|) trans- and gauche- molecules alternate 
and the intra-layer 0-distance (2.73 A) nearly equals the 
inter-layer 0-distance (2.72 A), so that the splitting is 
suppressed^^. Therefore, the observed OH-band splitting 
for bulk C16H33OH shows the presence of the 7-form. Ei- 
ther the whole bulk material exhibits the 7-form or there 
is a mixture of 7- and /^-crystallites. The latter case is fre- 
quently observed^^'^^: in fact, in the range of wavenum- 
bers from 1150 cm~^ to 950 cm~^, where C-C stretching 
vibrations are visible, we see indications for a superpo- 
sition of both forms (not shown). On the other hand, 
pore confined C16H33OH shows no OH-band-splitting at 
low temperatures. This reflects that the molecular ar- 
rangement doesn't transform in the monoclinic 7-form 
but remains in an orthorhombic structure, i. e. only the 
P-foim is present (compare Figs. [2] and [3|). This result is 
in agreement with the x-ray data presented above. Upon 
cooling, both the bulk and the confined hexade- 
canol pass from an hexagonal Ru- phase into a 
crystalline phase. The bulk material undergoes 
a stronger structural change, i. e. there is a mix- 
ture of the orthorombic /3- and the monoclinic 
7-form. The latter one is suppressed under con- 
finement, so that only the /5-form remains, which 



is quite similar to the hexagonal structure of the 
R//-phase: the fact that the crystahites have to fit into 
nanopores of irregular shape might favor the geometri- 
cally more simple /3-form^'^^(see Fig. [6j). 
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FIG. 8: (a) IR spectrum in the OH - stretching range for bulk 
C16H33OH. At lower temperatures the peak shifts to lower 
wavenumbers and then splits into two peaks, (b) Spectrum 
for confined C16H33OH, where no splitting is visible. 



Now let us turn towards the scissor-vibration of the 
CH2 groups (bending mode) that will give us informa- 
tion about inner-molecular and inter-molecular force con- 
stants. The spectra are shown in Fig. [TOl At first, we 
want to discuss the bulk material. At high temperatures 
(liquid state) a superposition of two peaks at 1467 cm~^ 
and 1460 cm~^ is observed. In the intermediate temper- 
ature range {Ru phase; see Fig. |4]) the intensity of the 
peak labeld " 1" increases strongly. At low temperatures 
(C phase; see Fig. [2]) this band splits up into two peaks. 
The latter transition can be clearly seen in Figs. [TTh ) 
and [T2b) . where we display the peak positions and inten- 
sities as a function of temperature. The results are sim- 
ilar to those obtained for the bulk state of n-paraffines, 
that apart from the missing OH-group are similar in their 
structure, i. e. that have the same CH2-backbone^^. In 
IR-spectra only one CH2-scissoring-band is observed at 
high temperatures, i. e. intra-molecular interactions of 
the CH2 -groups are too small to lead to a series of distinct 
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FIG. 9: (a) Wavenumber co/27vc of the OH - stretching peak 
vs temperature for bulk C16H33OH [compare with Fig. ^)]. 
Three different phases are visible: a) above 321 K, b) from 310 
to 321 K, where the peak position appears at lower wavenum- 
bers and c) below 310 K where the peak splits up into two 
peaks, (b) Wavenumber of the OH - stretching peak for con- 
fined C16H33OH (compare with Fig. [Sj)). The transition be- 
tween the C and the R// phase seems to be smeared in a 
range around T = 291 ± 5 K. The R/j-liquid transition does 
not affect the OH-stretching. 



peaks. The band splitting at low temperatures has been 
attributed to inter-molecular interactions (see Ref4^ and 
text below). 

Qualitatively, a behavior similar to that of the bulk 
state is observed for confined C16H33OH (see Fig. [TOb). 
In the high-temperature liquid phase two overlapping 
peaks are visible. The stronger one, i. e. that at higher 
wavenumbers, undergoes an increase in intensity at 
about 304 K (see Fig. [T2b). indicating the transition 
from the liquid phase to the Ru phase, while the 
secondary peak at lower wavenumbers gets weaker and 
finally disappears. At the second transition temperature 
of T = 291 K the remaining strong peak splits (see also 
Fig. [TTb). The separation is not as distinct as for bulk 
material. These transition temperatures, T = 304 K 
and T = 291 K (see Figs. Hb andUlD), agree weh with 
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those obtained via x-ray measurements (compare with 
Table ini). 
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FIG. 10: (a) IR spectrum showing the CH2 scissor-vibration 
for bulk C16H33OH at various temperatures, (b) IR spectrum 
showing the CH2 scissor-vibration of C16H33OH confined in 
mesoporous Si at various temperatures. 
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FIG. 11: Wavenumber of the CH - scissor peak vs temperature 
for (a) bulk C16H33OH and (b) confined C16H33OH (the peak 
labels refer to Fig.fTOl V 



In the following we want to analyze the dynamics of the 
CH2-groups in order to check whether it is affected 
by the geometric confinement, e. g. by an interac- 
tion with the pore surfaces, by the hmited num- 
ber of neighboring molecules (finite-size-effects) 
or by structural changes. In a first approximation we 
can assume that it is not affected by the stretching of 
the OH - groups. On the one hand, there is the scissor 
vibration, where the angle a between the two CH-bonds 
oscillates around its equilibrium value a = 109.47° (see 
Fig.[T3|). In addition, symmetric and asymmetric stretch- 
ing vibrations of the CH-bonds are observable (for the 
values see Table Hill) . Let fa and fd denote the respec- 
tive force constants. These can be calculated from the 
measured vibration frequencies using Eqs. ([26l) - ([28l) (see 
Appendix A; the difference in calculating fa via Eq. ([27|) 
or Eq. ([28]) is below 3.5% confirming that the inner- 
molecular coupling terms can be neglected). Table IIIII 
shows the results for the liquid and the Ru phase. Nei- 
ther the phase transition liquid R// nor geometrical 



confinement does markely affect the innermolecular con- 
stants. 

In the Ru phase the molecules rotate about their 
long axis, so that the primitive cell consists of only one 
molecule per layer (see Fig. |4]). Therefore, no splitting is 
observed. But in the C phase (below 310 K for bulk 
and below 291 K for confined C16H33OH), where the 
molecules are arranged in a herringbone structure, there 
are two molecules per layer in the primitive cell. So 
the symmetry of the arrangement allows a splitting of 
the scissoring band and obviously the molecular inter- 
actions are sufficiently strong that we are able to ob- 
serve a double peak (see above. Figs. [TOl and [TT]) . The 
strength of interaction depends on the distances between 
neighboring H- atoms of adjacent chains and can be ana- 
lyzed using a formalism developed by Snyder (see Ref4^ 
and Appendix B). In Fig. [Ml we have sketched the or- 
thorhombic lattice of the crystalline C16H33OH subcell 
(a view on the a-b-plane perpendicular to the molecules 
axis). In what follows we restrict ourselves to this /3- 
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bulk 


confinement 


Hquid 


Rii 


liquid 


Rii 


scissor [cm~-^] 


1467 


1467 


1467 


1467 


sym. stretch [cm~-^] 


2854 


2851 


2854 


2851 


assym. stretch [cm~-^] 


2927 


2921 


2924 


2918 


fd [N/m] 


455 


453 


454 


452 


U [N/m] 


56 ± 1 


56 ± 1 


57 =b 1 


57 ± 1 



TABLE III: Wavenumbers V — ijj j (27rc) (with a; being the an- 
gular frequency and c the speed of light) and resulting stretch- 
ing and bending force constants in the liquid and Rn phase 
of bulk and confined C16H33OH. has been evaluated using 
both the Eq. {|27| and the Eq. JSH)). The difference yields the 
specified uncertainty. are in units of N Ivn (see Eq. (|5]) in 
Appendix A and Ref.^^). To get in units Nm/radP one 
has to multiply fa with d^, where d — 1,09 • 10~^° m is the 
CH bond length. 




FIG. 12: Integrated intensity of the CH - scissor peak vs tem- 
perature for (a) bulk C16H33OH and (b) confined C16H33OH 
(the peak labels refer to Fig. [TOl ). 

form, that is characteristic for confined C16H33OH (a 
quantitative analysis of bulk C16H33OH is difficult due 
to the superposition of 13- and 7-form). Assuming that 
the inner force constant does not change at the phase 
transition, the intermolecular force constants /3^j can be 
evaluated from the observed splitting of the scissor band 
as described in Appendix B (see Eq. ([30]) ). The val- 
ues needed are the lattice parameters (see Table H]) and 
the herringbone angle between the projection of the 
backbone and the a-axis (see Fig. [14]). The latter one is 
determined via Eq. ([29]) and the measured intensities of 
the two CH2-scissoring-peaks. For confined C16H33OH 
we have la = 13.29 and = 22.23 yielding an angle of 
C = 37.7° (see Fig. ^ ioi T = 245 K). We display the 
intermolecular force constants in Table HVl For compar- 
ison, we also list literature values for an alkane, C23H48 
at 90 K, which have been evaluated in the same way^^. 
This alkane and C16H33OH exhibit a similar structure: 
The backbones of the molecules consist of the same CH2- 
units and both take the /3-form at low-temperatures. 
In addition, also the values of the lattice constants for 



w w 

H H 

FIG. 13: CH2 molecules with C-H bondlength d, H-C-H angle 
a and the resulting inner force constants fd and 

C23H48, a = 7.45 A and 6 = 4.96 A, are close to those of 
C16H33OH (see Tab.|I]). Due to this structural similarity 
the intermolecular distances listed in Table HVl are simi- 
lar, however, the respective force constants differ slightly 
by 10 to 20%. This is mainly due to the orientation of 
CH2-groups (the projection of the backbones on the a- 
b-plane) characterized by the herringbone angle (. For 
C16H33OH C = 37.7° holds, for the alkane C = 42°. 

This difference is probably due to the presence of po- 
lar OH-groups in C16H33OH that are strongly interacting 
and thus have an impact on the molecular orientation. 
The above comparison confirms once again that confined 
C16H33OH takes the /3-form in contrast to the bulk ma- 
terial (7- and /5-form). 

In order to assess the validity of our analysis, we 




FIG. 14: Lattice of C16H33OH in the C phase of the /3-form, i. e. a view on the a-b-plane perpendicular to the molecules axis 
(compare with Fig. [3]). Solid circles are H-atoms in the same plane, dashed circles are H-atoms lying in a plane above. 



also calculate the theoretical band splitting of the CH2- 
scissoring vibration and compare it with the measured 
values. Using the values from Table IIVI as well as 
Eqs. (|33|) and we get a theoretical value of AVcaic = 
8.1 cm-i for confined C16H33OH at T = 245 K. The mea- 
sured band splitting is AVmeas = 7.8 cm~^. Therefore, 
the experimental data is in good agreement with the the- 
ory. 



H-Atoms 


force 
constant 


confined C16H33OH 
(this work) 
C = 37.7° 


bulk C23H48 
(Ref^O) 
C = 42° 






distances (A) 


10~2^ Nm 


distances (A) 


IQ-^^ Nm 


3-2' 


/3,2' 


2.85 


-4.674 


2.79 


-5.837 


3-5 


/3,5 


2.59 


-2.775 


2.70 


-2.376 


3-6 


/3,6 


2.96 


-3.508 


2.94 


-3.920 


3-6' 


f3,6' 


2.96 


-3.508 


2.94 


-3.920 



TABLE IV: Distances between two neighboring hydrogen 
atoms and the resulting interaction force constants in the C 
phase confined C16H33OH (/3-form). ( denotes the herring- 
bone angle between the projection of the backbone and the 
a-axis (see Fig. I14|) . For comparison, we also list literature 
values for an alkane, C23H48 at 90 K, that also takes the /3- 
form^*^. The CH2-backbones of both molecules are similar, 
but C16H33OH exhibits an additional polar OH-group that 
affects the orientation angle 



IV. SUMMARY 

We have studied the structure and molecular dynam- 
ics of n-hexadecanol confined in nanochannels of meso- 
porous silicon and of bulk n-hexadecanol in their respec- 
tive phases (in the order of decreasing temperature: liq- 
uid, rotator R// and C). For this purpose we have per- 
formed x-ray and infrared- measurements. 

The transition-temperatures for confined C16H33OH 
are lower than for bulk C16H33OH (AT 20 K, see Ta- 
ble [III), addition, under confinement the phase tran- 
sitions are smeared, probably due to a distribution of 
pore diameters. Geometrical confinement does not affect 
the innermolecular force constants of the CH2-scissoring 
vibration (see Table Hill) but has an impact on the molec- 
ular arrangement. The R// phase of both bulk and con- 
fined hexadecanol is characterized by an orthorhombic 
subcell, where the chain axis are parallel to the layer 
normal (see Fig. (31). However, in the low-temperature C 
phase there is a fundamental structural difference. While 
bulk C16H33OH exhibits a polycrystalline mixture of f3- 
and 7-forms (see Figs. [2] and [3]), geometrical confinement 
favors a phase closely related to the /3-form: only crys- 
tallites with an orthorhombic subcell are formed, where 
the chain axes are parallel to the bi-layer normal. How- 
ever, the 7-form having a monoclinic subcell, in which 
the chain axis are tilted with respect to the layer normal, 
is suppressed. A reason for this might be the irregular 
shape of the nanochannels, into which the crystallites 
have to fit, favoring the formation of the geometrically 
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more simple and less bulky form^^^^ (see Fig. [6]). Since 
only the pure /3-form is present under confinement, we 
were able to evaluate the inter-molecular force constants 
of the CH2-scissor vibration. Also the orientation of the 



/3-crystallites has been determined: the molecules are ar- 
ranged with their long axis perpendicular to the pore 
axis. 
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Appendix A 

In this section we show how the innermolecular force 
constants of the CH2-groups can be evaluated using 
three characteristic vibration frequencies, that are easily 
measured: the scissor vibration as well as the symmetric 
and asymmetric CH-bond stretching. For this purpose 
we apply the Wilson FG - matrix method^^. We use the 
notation of Meister and Cleveland for the similar H2O 
molecule"^^ and perform the calculations in the same way. 

Fig. [13] shows a single CH2 - molecule. In the following 
we will neglect the influence of the neighboring molecules 
on this one. d = 1.09 A is the length of the C - H 
bond and a = 109.47° the angle between the two C - H 
bonds^^. This kind of molecule belongs to the €2^ point 
group. This means, there are two vibrations of type Ai 
(symmetric stretching and bending vibration) and one 
vibration of type B2 (asymmetric stretching vibration). 
The internal coordinates of this molecule are A(ii, Ad2 
and Aa. Adi and Ad2 mean changes in the bond length 
of the two C - H - bonds and Aa changes in the angle 
between the two bonds. Therefore we get three symmetry 
coordinates, two for Ai and one for B2. If we assume d 
being the equilibrium C - H bond length, then we obtain 
for the three symmetry coordinates: 




and are the transposes of r and R. With Eqs. ([7j) 

+ m 

r'fr = R'FR (9) 

The i?i's are linear combinations of the n's 

k (10) 
R = Ur 

Since the i^^'s are orthogonal and normalized, then 
U"^ = and 

r = U'R (11) 

= (U'R)' = R'U (12) 

This means with Eqn. (10) 

R'(UfUOR = R'FR (13) 
F = UfU' (14) 

The F matrix is 





Adi 


Ad2 


Aa 


Adi 


fd 


fdd 


dfda 


Ad2 


fdd 


fd 


dfda 


Aa 


dfda 


dfda 


d^fa 



The U matrix for type Ai is 



i?2 = Aa • d 
r. . 11 



i?3 = V^Adi- 



Ado 



(2) 
(3) 



Now, we have to calculate the F matrix, related to the 
potential energy, and the G matrix related to the kinetic 
energy. The potential energy can be written as 



2V = ^fikTiTk 

and with the internal coordinates 

2V =fa [(Adi)' + (Ada)'" 

+ faidAaf + 2/d„ (Adi + Ada) (dAa) 
+ 2fM (Adi) (Ada) 

Now we set di = d2 = d and write Eq. ([5]) as 

2V = J2f,iRjRi 
with Fji = Fij. In matrix form, Eqs. ^ and ([5]) become 



(4) 



(5) 



(6) 



and 



2V = r fr 
2V = R'FR 



(7) 
(8) 



and for B2 





Adi 


Ada 


Aa 


Ri 




^- 





R2 








1 


B2 


Adi 


Ada 


Aa 


R3 


^- 








So, for the type Ai the F matrix is 
F., = UfU' 





Fu ^ 




\ F21 


F22 J 





fd + fdd V2dfa 

V2dfa d^fa 



and for the B2 type 

(^33) 



{fd - fdd) 



(15) 
(16) 



The exact derivation of the G matrix shouldn't be shown 
here. It can be gleaned by Meister and Cleveland^^ . Only 
the most important steps shall be explained here. 
If only non-degenerate vibrations are present, the ele- 
ments of the kinetic energy matrix can be written as 



I 



(17) 
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where j and 1 refer to symmetry coordinates used in deter- 
mining the S vector, p refer to a set of equivalent atoms, 
a typical one of the set being t. fip is the reciprocal of 
the mass of the typical atom tp and gp is the number of 
equivalent atoms in the pth set. The S vector is given by 



lit) _ 



E 

k 



UjkSkt 



(18) 



where j, Ujk and have the same meaning as above. 
Skt can be expressed in terms of unit vectors along the 
chemical bonds and depends on the changes in the bond 
length or the angle between the bonds. So, the G matrix 
for the Ai vibration type has the form 



Gil Gi2 \ 
G21 G22 J 

IJ^H + Mc(l + COS a) 

lic^/l'^) sin a 



(19) 



+/ic(l-cosQ;) 
d 



and for the B2 vibration type 

GbsIG'ss) = {l^H + Mc(l - COS a)) 



(20) 



To determine the frequencies, one has to solve the equa- 
tion 



IGF- AEI =0 



(21) 



where \ = uo^ = {u27tc)'^ denotes the square of the angu- 
lar frequency. For the Ai type one gets the equation 



and for the B2 type 



A(FnGii 4 


- 2F12G12 - 


^11 F12 




Gil G12 


F21 F22 




G21 G22 



A.s — F^i^iG 



33 '-^as 







(22) 



(23) 



Eq. ([22|) can be separated with the Vieta expression^. 
Inserting the terms for the Fij and G^j, we obtain 

Ai +A2 = (/(i + /dd)[/ic(l + cosa) 

+ 2/c,[/ic(l - cos a) + fiH] - ^fda/^c sin a (24) 
Ai • A2 = [{fd + fdd)fa - 2/l]2/i^(2/ic + m) 
For Eq. ([23|) one obtains 

A3 = {fd -/(id)[Mc(l - cos a) (25) 

Neglecting the coupling constants fdd and fda allows 
to evaluate the innermolecular force constants using the 
measured wave numbers, Vd^sym = a/Ai /(^ttc), Va = 

\f}v2l (27rc) and Vd.asym — 

^f%l(2i^c). Then Eq. (f25l) 

yields 



772 

d^asym 



/J.c{l - COS a) + flH 



(26) 



Inserting this result into Eq. ([24|) yields 



fa = {2nc)^ 



• ^Ic l^c{'^ - cosa) + jiH 
2/i//(2/ic + Mi/) 



d,asym 



(27) 



There is a second possibility to evaluate i. e. by in- 
serting Eq. ([26]) into Eq. This yields 



fa = (2^c) 



—2 I —2 

2 ^ d^sym ^ a 



-2 _ /^^(l+cos a)-\-fiH 

d,asym ' /ic (1 -cos Q!)+/i// 



2(/ic(l - COS a) + jj^h) 



(28) 

Taking the measured wavenumbers listed in Table IIIII 
and the average angle between the CH-bonds, a = 
109.4°, as well as the masses of the atoms, 1/fic = 12ii 
and 1//J.H = lu {u = 1.6606 • 10-^^ kg), Eqs. ([26])-([28]) 
yield the force constants listed in Table IIIII The differ- 
ence in calculating fa via Eq. ([27j) or Eq. ([28|) is below 
3.1% confirming that the inner- molecular coupling terms 
can be neglected. 



Appendix B 

What follows is a summary of Snyder's derivation 
of the intermolecular force constants between the CH2 
groups of neighboring molecules that gives rise to a split- 
ting of the scissor band at low temperatures^. We show 
how this formalism can be applied to C16H33OH. An al- 
ternative description can be found in Ref.^^. 

In Fig.[T4lwe display a hexagonal subcell of C16H33OH. 
While Stein^ has taken only one pair of neighboring CII2 
into account to calculate the splitting of rocking and scis- 
soring bands, Snyder has shown that more pairs have to 
be included. When we consider the distances of H-atoms 
from the H-atom no. 3 (see Fig. [T4|) . then all atoms ex- 
cept no. 2', 6^ 5 and 6 have distances larger then 3.7 A. 
The internal coordinates ai are always half of the angle 
between the C - H bonds of a CH2 molecule. Solid circles 
are H-atoms in the same plane, dashed circles H-atoms 
in a plane above or below. 




FIG. 15: Lateral view at the long axis of the C16H33OH chain, 
X is the projection of the C - C distance on the c-axis of the 
crystal lattice 

Now, we want to write the positions of these five H-atoms 
as a vector. Fig. [15] shows the lateral view of a part of 
the C16H33OH chain. With values from Abrahamsson^ 
for Icc = 1.545 A and accc = 110.4°, we can calculate 
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the distance of the a-b-plane to the corresponding plane 
above or below with 

x = /ccsin(^^) = 1.2687 A 

Assuming that the hydrogen in the central plane has the 
c component 0, the hydrogen in the plane above has the 
component c = 1.2687 A. 



c 



a-b :plane 




plane i - 1 



FIG. 16: Lateral view at the long axis of the C16H33OH chain, 
d is the projection of half of the C - C distance on the a-b- 
plane of the subcell. 

The projection of the C - C bond in the a-b-plane is 
according to Fig. [16] 

d=^|icos(^)= 0.4409 A 

Taking the point for the lower left edge of the ab - plane 
the five atoms have the coordinates: 



Ho 



Ha 



H(\' = 



/ I — (icos(C) — I cos(a3 — Q 
I + (isin(C) — / sin(Qf3 — Q 

/ -dcos(C) +/sin(a2/ - f + 
h — (isin(C) — / cos(a2/ — f + C) 

V 

^ (icos(C) — I sin(a5 — f + C) 
(isin(C) + I cos(a5 - f + C) 
\ 1.2674 

^ (icos(C) + I cos(a6 — C) 
(isin(C) — /sin(a6 — C) 
^ 1.2674 

^ (icos(C) + / cos(a6/ — C) 

h + (isin(C) — I cos(a6/ — C) 
^ 1.2674 



with a and h being the lattice constants of the crystalline 
phase. We get the distances between the atom 3 and the 
other ones (see Fig. [M]) with 



ra,- = \H^ - Hj 



(28) 



where j = 2', 6^ 5, 6 is. The herringbone angle ( (see the 
upper left corner of Fig. [T4j) can be determined with the 
relation^ 



tan^C 



(29) 



where la is the integrated intensity of the scissoring 
mode, which is polarized in the a direction (higher 
mode at 1473 cm~^) and the one of the mode, which 
is polarized in the h direction (lower mode at 1462 cm~^). 

With the distances of two hydrogen atoms and Hj 
(=^3,j) [in our case 3 denotes the central H-atom (see 
Fig. [T4j) and j = 2', 5, 6, 6' the neighboring H- atoms that 
interact] we obtain the intermolecular force constants 

hj' 



dasdaj 



Vhh \ f dr\ rdr_ 



(30) 



where Vhh is the hydrogen repulsion potential intro- 
duced by Dows^^: 



Vhh = 1.2- 10"'^e 



10^-3. 52r 



(31) 



with r in A. 

The values of {O'^Vhh /dr'^)rij are obtained from 
Eq. 



Q^2 



1.486848 • lO-^e"^-^^'' 



(32) 



with /? in ^ = 10'^^^ = 1013 iv_ 

The measured intensity ratio (Eq. ([29|) ) allows us to 
calculate the distances Vij as well as the partial deriva- 
tives dr/dai (see Eq. (|28|) and above). Finally, by know- 
ing the intermolecular force constants fsj from Eq. (|32|) 
we can evaluate the band splitting of the scissoring vi- 
bration^^. For the angular frequencies 



T72 r-2 



f-L 



• Gf • {2/3,2' - 2(/3,6 + /3,60 - 4/3,5} 



holds ^ with 



(33) 



^ In Snyder's general theory the force constants for the scissoring 
vibration are denominated as (= /g^sO. fb (= h,6 + /3,60 
and /3 {= /3,5). 
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= ^Qlfic + QlfiH (34) 

Here l//ic = 12i^ and 1/ jiin = denote the masses of 
the atoms {u = LGGOG-IO-^^ kg), 1/Qr = 1.545-10-^0 m 
the C-C distance and 1/Qr = 1.09 • 10"^° m the C-H 
distance, so that Gf = (0.88825/ix) A-^ = 5.349-10'^^ (N 
m)-^s-^. For the band sphtting of the wave numbers we 
thus obtain 



Inserting the values of Table HVl i. e. (2/3^2' — 2 (/a, 6 + 
/3,60-4/3,5) = 15.784-10-2^ Nm for confined C16H336H, 
as well as the respective wave numbers, that we take from 
Fig. [m at low temperatures (labels la and lb: Vi-\-172 — 
2 • 146700 1/m) we get AV = 810 m = 8.1 cm"^ for 
confined C16H33OH. 



^2 



\2TrcJ 



Gf ■ {2/3,2' - 2( 



13,6 



/3,6' 



4/3,5} 



1^1 + V2 



(35) 



